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Starting with a model Hamiltonian, we study using the uniform expansion method conformational 
behavior of polyelectrolytes in the presence and absence of salt. The uniform expansion method 
yields all the important local length scales in the polyelectrolyte: the electrostatic blob size at large 
fraction of charges, the thermal blob size at low fraction of charges and the sizes of pearls (beads) 
and necklaces (strings) at intermediate fraction of charges. In the presence of salt, the electrostatic 
blob size depends on the ionic strength and increases with the increase in the salt concentration. 
We determine the salt concentration at which the pearl necklace intermediates dissolve and the 
nature of the transition changes from discontinuous to continuous. This critical salt concentration 
corresponds to the length scale where the Debye screening length is of the order of the necklace 
length. 

PACS numbers: 36.20.-r, 82.35.Rs, 64.60.-i 



I. INTRODUCTION 



Polyelectrolytes are charged polymer with ionizable groups that dissociate ions in solution, leaving the charges of 
one sign bound to the chain and counterions in solution 1]. Inspite of the hydrophobic backbone, the presence of 
the charged side groups make these molecules water soluble, a property which is very important for biological and 
technological applications % l3( . The competition between the short range attractive interactions due to poor solvent 
conditions and the long range repulsion between like charges gives rise to structures not seen in neutral polymers. The 
rich conformational behaviour of polyelectrolytes can be finely tuned by changes in the ionic strength, pH, temperature 
and condensation agents J|. However, the presence of the long-range electrostatic interactions in combination with 
the polymer elasticity make these system difficult to handle both analytically and computationally. 

Simulations on polyelectrolytes in poor solvents show a first order transition from a globular state at low fraction 
of charges to an extended state at large fraction of charges. The transition proceeds via the formation of metastable 
pearl-necklace intermediates - locally collapsed globules joined by narrow strings. Pearl-necklace structures have been 



seen in simulations both in the absence 



a, M, and presence of counter-ion condensation [7|, ||J |9|, llfj . The splitting 



of a charged globule to form a pearl-necklace structure is similar to that of the Rayleigh instability of a charged 
liquid droplet, where system gains energy by splitting into smaller dro plet s placed far apart 
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121 ]. Scaling theories 



16[ . However, there are very few 



have provided a deeper insight into the structure of pearl necklaces 
analytical theories, mainly based on the Feynman variational method, that confirm the existence of pearl-necklaces 
in the absence of salt and counter-ion condensation 
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181 ]. One reason for the limited literature is that analytical 



models require explicit contributions from chain entropy, short range two-body attractive, three-body repulsive and 
long-range screened Coulombic interactions. The inclusion of all these terms is a formidable task in a theoretical 
description. However, there are methods like the uniform expansion method, which provide systematic and useful 
way of looking at the systems with many different length scales 
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In an earlier shorter communication, one of us used the uniform expansion method in the absence of salt to establish 

rn 

that the globule-to-rod transition is brought about via intermediate pearl-necklaces, which are metastable 22]. The 
transition was found to be first order in nature. The critical fraction of charges where the first-order transition 
occured was estimated using Maxwell equal area construction. Although the discontinuous nature of transition was 
established, the mean size of pearl-necklaces and several other important length scales in polyelectrolytes could not be 
calculated. Also, how the presence of salt can change the nature of transition and dissolve pearl-necklace structures 
was not considered. 

In this work, we present a detailed calculation for polyelectrolytes in poor solvents in the presence and absence 
of salt. We start with a model Hamiltonian for polyelectrolytes in poor solvents incorporating entropic, screened 
electrostatic, two body attractive and three body repulsive terms. We first obtain the mean size of the electrostatic 
blobs, the thermal blobs, the beads, the strings and pearl necklaces using the uniform expansion method in the 
absence of salt. The present study exploits the perturbative nature of the uniform expansion method to estimate the 
important length scales of a polyelectrolyte in a poor solvent. The Hamiltonian based approach adopted here allows 
for a complete investigation of the phase transition between different conformational states. The details are presented 
in Section III. In the presence of salt, the screened electrostatic interactions are characterized by the Debye screening 
length. Depending on the ratio of the Debye screening length and the size of the polymer, the nature of rod-to-globular 
transition changes from discontinuous to continuous. In Section IV, we study in detail the nature of transition as a 
function of the salt concentration. We determine the salt concentration which dissolves the pearl-necklace structures 
and brings about a direct transition. 
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II. MODEL 

The Hamiltonian of an infinitely dilute solution of polyelectrolytes in a poor solvent in the presence of salt is given 

by 

H = Ho + H2+H 3 + H c , (1) 



where 



H= 2VJ d *Hw ~2i ds J Q ds ' 5 lr e / ds J o ds 'j Q ds"8[r(s)-r(s')}8[r(s')-r(s")] 

where the polymer conformation is decribed by the radius vector r(s) at distance s from one end of the chain; b is 
the Kuhn length and N is the number of monomers. Tig describes the entropic elasticity of the chain; H.2 1 T~i3 and 
Ti. c account for the two-body attractive, three-body repulsive and the screened coloumbic interactions respectively; 
V, w and v c represent the strengths of the two-body, three-body and electrostatic interactions respectively. The 
temperature dependence of v is given by v — r6 3 , where t = \T — 9\/9 is the reduced temperature and 9 is the 
theta temperature [22j |. The sign of the two-body interaction term is positive under good solvent conditions and 
negative under poor solvent conditions, k is the inverse screening length which depends on the ionic strength by 
r^ 2 = k 2 = AttIb Yli c i1i i where rjj is the Debye screening length; Cj is the concentration of the ith ion with charge 
qi] Ib — e 2 /eksT is the Bjerrum length defined as the distance at which the electrostatic interaction between the two 
elementary charges is equal to the thermal energy fcgT and e is the dielectric constant of the medium. The strength 
of the electrostatic interaction is given by v c = f 2 ls, where / is the fraction of charges along the chain backbone. It 
is to be noted that the fraction of charges are assumed to be too small to bring about counterion condensation. 
In terms of the Fourier transform, the interaction energy can be written as 

r N r N p >3 k »k-|r(s)-r(s')| 



The calculation of the mean size of the polyelectrolyte using the above Hamiltonian is a difficult task. A self-consistent 
perturbative approach of Edwards and Singh known as the uniform expansion method can be used to calculate the 



polymer size 



approa 

mm 



21] . The uniform expansion method is based on defining a new step length, b\ 3> b, such that the 
mean square end-to-end distance of the chain in presence of the excluded volume remains Gaussian and is governed 
by (R. 2 ) = Nb\. This can only happen if the original Hamiltonian Ho = ^f? dsr(s) 2 is replaced with the reference 
Hamiltonian Hi = J Q dsr(s) 2 . This step can be carried out by adding and subtracting Hi from Ho to give the 
following expression for the mean square end-to-end distance: 



<R 2 ) = 



J V[r{s)}R 2 expj-THi + {Hp - Hi) + H 2 +H 3 + H c )} 
J V[r(s)} exp[-(Wi + (H Q - Hi) + H 2 + H 3 + H c )\ ' 



(6) 



One can now expand (R. 2 ) in a perturbative series about the reference Hamiltonian Hi such that to the first order 
correction in v one obtains the following variational equation for the unknown parameter bi: 



<R 2 ) 1 (Ho -Hi+H 2 + H 3 + H c ) 1 - (H 2 (H Q -Hi+H 2 + H 3 + H c )) 1 = 0, 



(7) 



where the average is with respect to the following probability distribution: 



"P[r(s)] a exp 



3 f N 

^Jo dsHs)2 



Since the probability distribution is Gaussian, the evaluation of the averages in Eq. ((7]) is simple. The details of the 
calculation are presented in Appendix A; here we simply present the result: 



2 tNV* f 2 l B N*/ 2 F(P) 1 

1 — a 1 ; -rr-p. 1 T = U. 

^ 3 ba/3 5 ' 2 ™ 6 



O"' 



a" 



(8) 



In the above equation, a = R/Ro is the expansion (or compression) factor, which characterizes the change in the 
polymer size with respect to the unpertubed size Rq — N 1 / 2 ^ v c — J 2 Ib', v — rb 3 ; w — b 6 ; (3 — n 2 Nbf/6 = {R^) /rjj, 
where (3 represents the dimensionless ratio of the mean square polymer size with respect to the square of the Debye 
screening length ra- is given by 



F{fi) = e^l - 0(v^)](l - 2/3/3 + f3 2 /6) + " 0/3 - 1, 



(9) 



5 

where <j) is the probability integral. The first two terms in Eq. ([8]) account for the entropy of the polymer; the third 
term is due to the excluded volume interactions between chain segments which is attractive in nature; the fourth 
term accounts for the electrostatic repulsion; the last term represents the three-body interactions. For convenience, 
we have ignored all the numerical coefficients. 



III. POLYELECTROLYTE CHAIN IN A SALT-FREE SOLUTION 



The limit of /3 — > in Eq. © represents the case of a salt-free solution. In this limit F(/3)//3 5 / 2 
reduces to 



1 and Eq. 



2 tN^ fH B N*l* 1 
I- a ; 1 1 g = 0, 

a a ba or 



(10) 



The above equation was solved analytically in Ref. 



22 1 in two opposite limits. In a poor solvent at low fraction of 
charges, a < 1, tN 1 / 2 ^> 1 and / 2 is-/V 3 / 2 jb <C 1, the attractive interactions dominate, and the comparison of the 
third and the last term gives a globule of size R g i b ~ fer -1 / 3 ^ 1 / 3 . At large fraction of charges, q»1, tN 1 / 2 <C 1 and 
PIbN 3 / 2 jb 1, the electrostatic repulsion overcomes the attractive interaction. As a result, the comparison of the 
second and the fourth term gives a = (/ 2 ?B/6) 1/ ' 3 A r1 / 2 resulting in a rod-like state, R ro d ~ bN(f 2 ls/b) 1 ^ 3 . Although 
the approach used in Ref. 



22j gives the right prediction for the globular and rod-like state, it can not predict the size 
of pearl necklaces. This is because the formation of pearl necklaces involves contributions from entropic, steric and 
electrostatic interactions, and therefore requires comparison of more than two terms. 

Strictly speaking, the uniform expansion method is a perturbative ap pro ach, where all the interaction terms other 
than the entropic contribution should be considered as a perturbation 
more than two terms analytically. In what follows, we make use of the perturbative nature of the uniform expansion 
method to calculate the local length scales 2l|, which are then used to estimate the pearl necklace size. In terms of 
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2l|. Such an approach can compare 



the blob picture, the above scaling can be reconsidered in the following limiting cases. 



A. Polyelectrolytes in poor solvents at low fraction of charges 



At low fraction of charges, the electrostatic interaction term (the fourth term in Eg. (fit)]) ) can be ignored. In this 



limit, an important local length scale is the thermal blob size 



231 ] . For sizes less than the thermal blobs £x , monomers 



1/2 

do not feel the attractive interactions due to the poor solvent conditions and follow gaussian statistics, i.e., £t ~ N T 6, 
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where Nt is the number of monomers in the thermal blob. For sizes larger than £r, the thermal blobs are collapsed 
to form a compact sphere, i.e., R g i b ~ {N/NtY^^t- To get the number of monomers in the thermal blob, let us first 
consider a region close to the O temperature, where r ~ 0. In this region, Eq. (fTO)) suggests that a ~ 1, and the chain 
follows Gaussian statistics with size given by i?o ~ N 1 ' 2 b. As the solvent quality becomes poor, the effects of the 
attractive interactions are realized at long length scales. However, at short length scales the chain remains ideal with 
the size dictated by £r. This implies that the excluded volume term (the third term in Eq. (|10p ) is a perturbation 
to the gaussian size a ~ 1 as long as tN 1 / 2 <C 1. The latter gives the number of monomers Nt ~ 1/t 2 inside the 
thermal blob of size given by 

£r ~ N T /2 b ~ 6/t (11) 

For N < Nt, the statistics is gaussian and the size is dictated by £t- For N > Nt, the attractive nature of the 
interaction term dominate and the thermal blobs are space filling with the size given by R g i b ~ (htY^^t, where 
tit = N/Nt are the number of thermal blobs. The overall size scales as 

Rglob ~ bT- l / Z N X / Z , (12) 

which recovers the size of a polymer globule at low fraction of charges. 

B. Polyelectrolytes in poor solvent at large fraction of charges 

At large fraction of charges, the electrostatic interactions dominate over the attractive interactions, and the latter 
can be ignored. This limit corresponds to theta solvent, where electrostatic blobs are linearly arranged. In this limit, 
an important local length scale is the electrostatic blob size, £, e iec- For sizes less than the electrostatic blobs, £eieo 

1/2 

monomers do not feel the repulsive electrostatic interactions and follow gaussian statistics, i.e., £ e z ec ~ N e [ ec b, where 
N e i ec is the number of monomers in the electrostatic blobs. For sizes larger than £ e / ec , the electrostatic interactions 
dominate and blobs are linearly arranged to result in an extended state, i.e., R ro d ~ (N / N e i ec )£ e i ec . In this limit, the 
electrostatic interaction term (the fourth term in Eq. IjlOp) is a perturbation to the gaussian size a ~ 1 as long as 
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f 2 l B N 3 / 2 /b < 1. Thus the number of monomers N e i ec that retain the gaussian statistics are given by 

N e l ec ~ {f 2 l B /b)-^ 3 , (13) 



1 /2 

such that for N < N e i ec the statistics is ideal with the size £, e iec ~ NJ ec b given by 



Ce/ec ~ (flB/b)- 1 ^. (14) 

Therefore, a short chain follows gaussian statistics. For TV > N e i ec , the electrostatic blobs are arranged linearly to give 
the mean size R ro d = (N/N e i ec )^ e i ec , where n e / ec = N/N e i ec are the number of electrostatic blob at a given fraction 
of charges. Therefore, the overall size is given by 

Rrod ~ bN(f 2 l B /b)^ 3 , (15) 

which recovers the polyelectrolyte rod size at high fraction of charges. In the limit where / ~ (l B /b)^ 1 ^ 2 N~ 3 ^ 4 , the 
polyelectrolyte rod size R roc i crosses over to the gaussian size i?o- 

C. Polyelectrolytes in poor solvents at intermediate fraction of charges 

To understand the behavior in the intermediate regime, Eq. pop can be solved numerically. The numerical solution 
is shown in Fig. (1), which is a plot of the dimcnsionless size a as a function of the dimensionless electrostatic energy 
/ = f 2 l B N/b. The figure shows an intermediate regime where the curve folds back on itself. This region represents a 
bistable equilibrium between the collapsed globules and elongated strings that can coexist. This is the region where 
pearl- necklace structures are formed. In the intermediate region between f min and f max , there are three curves: the 
two curves increases with the increase in the fraction of charges and represent the metastable pearls and necklaces; 
one curve decreases with the increase in the fraction of charges and represent the unstable region. 

A detailed analysis of this plot was presented in Ref. 



22j , where the Maxwell equal-area construction was used to 



determine the critical fraction of charges where the abrupt transition (spontaneous Rayleigh splitting) occured. The 
pearl-necklace structures formed at the intermediate fraction of charges are the result of the balance between the long 
range repulsive electrostatic interactions and the short range attractive interactions. In terms of the above scaling, the 
intermediate regime corresponds to f — f 2 l B N/b ~ r. This gives the critical fraction of charges, f c ~ {rb/lsN) 1 / 2 , 



8 

where the abrupt transition is expected to take place. For r = 1.0 in Fig. (1), the scaling predicts the transition at 
f c = f^lsN/b ~ 1.0. This is very close to the value of f c = 0.9 obtained from the Maxwell equal area construction 
used in Fig. (1). Therefore, at the intermediate fraction of charges, the number of monomers inside a Rayleigh blob 
5| can be obtained by comparing the third and the fourth terms in Eq. (TIT))) . Since the statistics at the local length 
scale is gaussian, a ~ 1, the number of monomers inside the blob are given by 



N R ~Tb/f 2 l B . (16) 

A globule splits into smaller globules as soon as the electrostatic energy, fc^T f 2 lBN 2 /%/ 3 £t, overcomes the surface 
energy, ^Tn^ 3 5|. The balance of the two terms gives the number of monomers in a Rayleigh blob, Nr ~ rb/ f 2 l B - 
The latter expression is the same as Eq. (|16[) . This is because the surface energy considerations are already there 
implicitly in the present formalism in terms of the fundamental interactions. 

It is clear from Fig. (1) that in the intermediate regime between f min and f max , the size of a pearl-necklace 
structure is determined by the size of the pearls (beads dbead) and necklaces (strings l s t r ) that contribute to the 
overall structure. The size of the pearl- necklace structure can be calculated in two limits: when db ea d 3> Istri the 
pearl- necklace structure is dominated by beads and corresponds to the regime f min ~ / < f c , represented by a dashed 
line on the left hand side of f c ; the opposite limit of l s t r ^> dbead corresponds to the regime f c < f ~ j max - It is 
worth mentioning that the size of pearl- necklaces follow universal scaling in terms of N /Nr only close to the globular 
and rod-like states. This is because the globular and rod- like regions are separated by a first-oder transition [5(. 

In what follows, the size of pearl-necklaces is estimated in the limits mentioned above. With the increase in the 
fraction of charges, the collapsed globule with size given by R g iob ~ (N/Nt) 1 ' 3 ^ is expected to break down into 
small sized globules (beads) , the size of which is given by 

dbead ~ {Nr/NtY'^T ~ N R /3 T-^b ~ (b/fl B ) 1/3 b. (17) 



It is clear from Fig. (1) that for f m i n ~ / < f c , the size of the bead dbead is much larger than the length of the string 
I str , i.e., dbead ^ htr- ln this limit, for N > Nr the electrostatic interactions determine the overall size 



Rp - (N/N R )d b ead ~ (f 2 l B /b) 2/3 bN/T. 



(18) 



This gives the size of a pearl- necklace chain in the limit dbead 3> l s tr- The size obtained in this limit is similar to the 
elongated globule picture of Khokhlov 24j . It is to be noted that the number of pearls are given by n p ~ N/Nr ~ 



(f 2 ls I 'rb)N . By changing the fraction of charges / or r, a globule splits into different (i nteger) number of pearls 
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14j . a phase diagram 



n p , resulting in a cascade of abrupt transitions between necklace globules. In References 
between the fraction of charges / 2 ?b/& and r show a cascade of transitions with different number of pearls, which 
occur because of the change in either / or r . When a globule splits up into progressively large number of smaller 
globules (pearls), the total size of the pearls decreases while the total length of the strings increases. In the context 
of present work, the region between f min and f max represents a cascade of transitions, where a globule splits into 
different number of pearls sizes as soon as / > j ' min . The decreases in the overall size of the pearls is a reflection 
of the increase in its number. The intermediate regions in Figs. (1) and (3a) represent two such instances where a 
cascade of transitions occur because of the change in / and r respectively. 

In the limit f c < / ~ f max i the length of the string l str is much larger than the size of the string dbead, i-e., 
l s tr S> dbead- On the local length scale, the length of the string is dictated by the gaussian statistics given by 



l str ~ {Nr/Nt) 1 ' 2 ^ ~ {Tb/fl B y/ 2 b. 



(19) 



The above scaling can also be obtained be comparing the surface energy of the string, kBTl str /^T, with the electrostatic 
repulsion between the pearls, fegT 'f 2 lBN R /l s t r [5]. On large length scale, the electrostatic interactions dominate and 
the overall size is given by 



R n ~ (N/N R )l str ~ (flB/br^bN, 



(20) 



which yields the pearl- necklace size in the limit l s t r 3> dbead- 

The size of the bead, string, pearl-necklace and the number of pearls calculated from the uniform expansion method 
is in complete agreement with the previous scaling theories 



13 



14 



15 



16[. The advantage of the present theory is 



that one can easily include the effects of salt. The details arc presented in the next section. 
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IV. POLYELECTROLYTE CHAIN IN THE PRESENCE OF SALT 

The limit of (3 ^> 1 in Eq. §§§ represents the case of high salt concentration. In this limit F(f3) — » /3 3 / 2 and Eq. ([5]) 
reduces to 

1-a 2 —+ J " + — = 0, 21 

a 3 fea/3 a 6 v ; 

Since (3 = K 2 Nb 2 = K 2 a 2 Nb 2 , the above equation reduces to 

1 - a 2 — + \ 3 B 3 2 + — = 0, 22 

In a theta solvent, the electrostatic term is a perturbation (a ~ 1) as long as ^IbN 1 / 2 /b z n 2 <C 1. This gives the k 

1/2 

dependent electrostaic blob size £ e *,sait ~ salt b, where 

iV ei , sait ~ (& 3 k 2 // 2 Z b ) 2 , (23) 
is the number of monomers inside the electrostatic blob of size given by 

Usalt ~ b 4 K 2 /fl B . (24) 

In a theta solvent, the third term in Eq. (|2ip can be ignored and the fourth term represent the excluded- volume type in- 
teraction term. In this limit, the comparison of the second and the fourth terms yields R sa , S ait — (/ 2 /b/6 3 k 2 ) 1 / 5 ^V 3 / 5 6. 
The latter represents a self avoiding walk with an effective persistence length given by 

l P ,sait ~ (f 2 l B /b 2 K 2 ) 1/5 b. (25) 

With the increase in the salt concentration the K-dependent electrostatic blob size will increase until £ e l,salt becomes 
of the order of £ e / ec . The comparison of the two gives the minimum value of the salt concentration that influences the 
electrostatic blob size, 



Kmin ~ (fh/b) 1 ^- 1 . 



(26) 
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For k < K m in, the charges on the chain backbone are unscreened and the size is given by R ro d ~ bN^pls/b) 1 ^ 3 . In 
the opposite limit, k > K m in, some of the charges are screened and the size is given by R sa , S ait — {f 2 lB /b 3 k 2 ) 1 /^ N 3 / b b. 

The blob size increases with increase in salt concentration until K = Kg. At this salt concentration, the effective 
excluded volume, v* = r — f 2 lB/b 3 K 2 , becomes zero and 

kb ~ (f 2 l B /b 3 T)^ 2 . (27) 

At k — Kg, the electrostatic blob size £ e l,salt becomes equal to the thermal blob size £t = b/r and the chain behaves 
like an ideal chain. For k > Kg, the effective excluded volume becomes negative due to which the attractive interaction 
dominates and the thermal blobs of size £t collapse to form a globule of size Rgiob ~ br^^N 1 / 3 . Thus the pearl- 
necklace intermediates are dissolved at the salt concentration given by k = Kg. 

The comparison of Eq. (|27[) with Eq. (fTO|) shows that Kg ~ l/l s tr implying ro ~ l str . The latter shows that the salt 
concentration that dissolves the pearl necklace structure corresponds to the case where the Debye screening length is 
of the order of the necklace (string) size. 

To understand these results quantitatively, Eqs. ijHJ) and © can be solved numerically. The result is presented in 
Figure 2, which is a plot of the dimensionless ratio of the chain size a as a function of the dimcnsionless inverse sreening 
length p. The curve shows a slow decrease in the chain size with the increase in the salt concentration followed by an 
abrupt collapse to a globular state. As suggested by the above scaling arguments, the abrupt transition is expected 
to take place at salt concentration given by Kg ~ (f 2 ls /b 3 T) x ^ 2 . In terms of the dimcnsionless parameters (3, the 
transition is expected to take place when f3 c — KgNb 2 /6. For r = 0.5, plsjb — 0.1 and N = 100, the scaling predicts 
the transition at f3 c = 3.3. This is very close to the value of /3 C = 3.0 in Fig. (2) obtained by solving Eqs. © and © 
numerically. 

Figures 3a-3c are the plots of the dimcnsionless size a as a function of the reduced temperture r at three different salt 
concentrations characterized by Kb = 0.17,0.35,0.42 respectively, obtained by solving Eqs. © and numerically. 
For the values of N = 100, f 2 ls/b = 0.05 used in Fig. 3a-3c, the minimum salt concentration required to influence 
the polyelectrolyte size is given by K m i n b ~ (f 2 lB /b) 1 ^ 3 ~ 0.37. Fig. 3a corresponds to the case where k < K m i n . As 
a result, the charges on the polyelectrolytes are completely unscreened and the first-order transition proceeds via the 
formation of an intermediate pearl-necklace structure; Figure 3b represents the intermediate regime where k ~ K m j n 
and some of the charges on the chain are screened. In this limit the fraction of charges on the chain are too small 
to form the metastable pearl-necklace structure and the transition is weakly first order. As before, the transition is 
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expected to take place when k ~ Kg ~ (/ 2 ^b/t6 3 ) 1 / 2 . For f 2 ls/b = 0.05 and k& = 0.35, the scaling suggests that the 
transition takes place at r c ~ 0.4. This is very close to the value of t c = 0.37 in Fig. 3b; Figure 3c represents the case 
where k > K m i n and the charges on the chain backbone are completely screened. In this limit the transition becomes 
smooth and continuous second-order in nature. 

V. CONCLUSIONS 

We have used the uniform expansion method of Edwards and Singh to study conformational transitions of poly- 
electrolytes in poor solvents in the presence and absence of salt. The presence of the additional long-range repulsion 
due to the electrostatic interactions between charged monomers destabilizes the globular structure as the fraction of 
charged monomers along the chain backbone is increased and leads to the formation of pearl-necklaces at intermediate 
fraction of charges . The uniform expansion method allows us to determine the the size of the pearls (beads) , necklaces 
(strings) and pearl-necklace intermediate, and shows a first order transition between a globular state to an extended 
state via metastable pearl-necklaces. 

In the presence of salt, the electrostatic blob size depends on the salt concentration and is found to be larger 
than that of the unscreened case. We determine the minimum salt concentration, K m i n , at which the size of the 
electrostatic blob is influenced. For k < K m in, the size of the polyelectrolyte rod is governed by the unscreened 
electrostatic interactions. For k ~ Kmim 011 the other hand, some of the charges are screened and this results in the 
increase in the electrostatic blob size amounting to decrease in the polyelectrolyte size. The electrostatic blob size 
increases with the increase in the salt concentration until it becomes equal to the thermal blob size. We determine the 
salt concentration, Kg, at which the thermal blob size is equal to the electrostatic blob and the polyelectrolyte behaves 
like a neutral ideal polymer. This is the concentration at which the pearl-necklace intermediates are dissolved and the 
transition proceeds directly. For k > Kg, most of the charges on the chain backbone are screened and the attractive 
interactions dominates resulting in the chain collapse to form a globule. Our scaling predicts that the pearl-necklace 
intermediates dissolve when the Debye screening length is of the order of Z str . In a future publication, the present 
formalism will be extended to include the effects of counterfoil condensation [261 ] . 
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APPENDIX A: UNIFORM EXPANSION METHOD 



The detailed calculations of the terms containing the entropic term, (R- 2 ) 1 ("Ho — ~ (R 2 (7^o — an d the 



two-body interaction, (R 2 ) 1 (H2 + ^c))i — (R 2 ('^2 + H C )Y, are fairly standard and can be found in References 



and 



20]. 



19] 



{R 2 ) 1 (Ho-H 1 ) 1 -{Tl 2 (Ho-H 1 )) 1 
{R 2 ) 1 (H 2 ) 1 -{R 2 (H 2 )) 1 
<R 2 ) 1 (W C ) 1 -(R 2 (W C )) 1 



Nh 4 



1 1 

2vh 4 f°° f N 
— xo / ak as 
9(2^)2 J Q J Jtl 

9(2 



ds '( s _ s ')2 fc 4 e -fc%2( S - S ')/6 5 



A' 4 



-fe 2 6i 2 (s-s')/6 



(fc 2 + K 2 ) 



(Al) 
(A2) 
, (A3) 



The uniform expansion method of Edwards and Singh does not include the three body interaction term. A detailed 
calculation of the three body interaction term is given in Ref. 23. In what follows, we present a few important steps 
required to calculate the three body interaction term, (R 2 ^ (^3)1 — (R- 2 ('^3)) 1 - Let us first calculate (R 2 ?^)^ 



pN rN rN 

(R 2 H 3 ) 1 = - ds ds' ds" (R 2 6[r( s ) - r(s')]5[r(s') - r(s")}) 
Jo Jo Jo 



(A4) 



In terms of the wave vectors k and q, the above equation can be rewritten as 



"' 1 ^ f °° fN ds [ N ds' / ,JV ds "/ R a e ik-kW-r(-')l e iq-W.')-r(V')| 



(2tt)3 (2tt)3 J J 



(A5) 



To calculate the quantity of interest, it is important to expand the mean square end-to-end distance and the 
probability distribution in terms of the internal coordinates. Since the probability distribution is Gaussian, the 
averages can easily be calculated to give the following expression: 



(s - s') 2 kV ( s < - s ") 2 qV 



9 



e -k 2 (s-s')b 1 2 /e e - tl 2 (s'-s")b 1 2 /e 



(A6) 



One can repeat the above procedure to calculate (R 2 ) 1 (^3)1, which is given by 



/°° ^31, r°° fN rs rs 

7^4 / W% / ds ds' ds''Nb x *e-^-WlSe-^'-°''^/\ (A7) 
-00 ( Zn ) J -oo V Zn ) Jo Jo Jo 
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After substracting Eq. (A7) from Eq. (A6), one obtains 



(R 2 ) x (H 3 ) 1 - (R 2 H 3 ) 1 = 



wfl-^wfL ds L ds 'L ds " 



{s - s') 2 k 2 6i 4 (s' - s") 2 qV 



e -k 2 ( S - S ')ftl 2 /6 e -q 2 (s'-s")bl 2 /6 



(A8) 



The above equation can be rewritten as 



(R 2 ) 1 (W 3 ) 1 -<R 2 (W 3 )) 1 



4 />OG />oo 



N 



I dk I dq I ds I ds' I ds" [(s - s'fk^q 2 + (s' - s") 2 k 2 q 4 ] 
9(27r) Jo Jo Jo Jo Jo 

e -k 2 b 1 2 (s-s')/6 e -q 2 b 1 2 (s'- S ")/6 



(A9) 



The last integral diverges as q — > oo. The divergence can be removed by introducing an upper cut-off for the wave 
number q, the details of which are discussed in Ref (23). The integrations in Eqs. (A2), (A3) and (A9) can easily be 
carried out, the results when substituted in Eq. |(7j) gives the following variational equation: 



3 / 2 v l/2 



wN 



= 0, 



(A10) 



where F(j3) =e [l- </>(v^)](l - 2/3/3 + /3 2 /6) + - (3/3 - 1; <f> is the probability integral; (3 = Hi 2 Nb 2 /Q and 

k is the inverse screening length. The above equation can be written in a dimensionless form by dividing it by Nbi 2 
and defining a — b\/b. The final form of this variational equation is given by Eq. I|8|). 



[1] J.-L Barrat and J-F Joanny, Adv. Chem. Phys. 94, 1 (1996). 

[2] M. Hara, Polyelectrolytes; Marcel Dekker: New York, 1993. 

[3] K. S. Schmitz, Macroions in solution and colloidal suspension, 1st ed.; VCH Publishers: New York, 1993. 

[4] A. Kiriy, G. Gorodyska, S. Minko, W. Jaeger, P. Stepanek and M. Stamm, J. Am. Chem. Soc. 124, 13454 (2002). 

[5] A. V. Lyulin, B. Duenweg, O. V. Borisov and A. A. Darinskii, Macromolecules 32, 3264, (1999). 

[6] P. Chodanowski and S. Stoll, J. Chem. Phys. Ill, 6069 (1999). 

[7] U. Micka and K. Kremer, Europhys. Lett. 49, 189 (2000). 

[8] H. J. Limbach and C. Holm, J. Chem. Phys. 114, 9674 (2001). 

[9] C. Holm, H. J. Limbach and K. Kremer, J. Phys. Condens. Matter 15, S205 (2003). 

[10] H. Limbach, C. Holm and K. Kremer, Europhys. Lett. 60, 566 (2002). 



15 

[11] Y. Kantor and M. Kardar, Europhys. Lett. 27, 643 (1994). 

[12] Y. Kantor and M. Kardar, Phys. Rev. E 51, 1299 (1995). 

[13] A. V. Dobrynin, M. Rubinstein and S. P. Obukhov, Macro-molecules 29, 2974 (1996). 

[14] A. V. Dobrynin and M. Rubinstein, Prog. Polym. Sci 30, 1049 (2005). 

[15] H. Schiessel and P. Pincus, Macromolecules 31, 7953 (1998). 

[16] H. Schiessel, Macromolecules 32, 5673 (1999). 

[17] G. Migliorini, N. Lee, V. Rostiashvili and T. A. Vilgis, Eur. Phys. J. E 6, 259 (2001). 

[18] F. J. Solis and M. Olvera de la Cruz, Macromolecules 31, 5502 (1998). 

[19] M. Doi and S. F. Edwards, The Theory of Polymer Dynamics; Clarendon Press: Oxford, 1986. 

[20] S. F. Edwards, P. J. Singh, Chem. Soc. Faraday Trans. 2, 75, 1001 (1979). 

[21] P. G. Higgs and J. F. Joanny, J. Chem. Phys 94, 1543 (1990). 

[22] A. Dua and T. A. Vilgis Europhys. Lett. 71 49 (2005). 

[23] P. -G. de Gennes, Scaling Concepts in Polymer Physics; Cornell University Press: Ithaca, NY, 1985. 

[24] A. R. Khokhlov, Journal of Physics A 13, 979 (1980). 

[25] A. Dua and T. A. Vilgis Macromolecules 40, 6765 (2007). 

[26] M. Muthukumar, J. Chem. Phys. 120, 9343 (2004). 



16 



Figure Captions 



Figure 1. Variation of the dimcnsionlcss chain size a as a function of the fraction of charges / = f 2 lBN/b for 
t = 1.0 and TV = 100. The region between f min and f max represents a bistable equilibrium between pearls and 
necklaces. The Maxwell equal area construction has been used to locate f c , the critical fraction of charges where 
spontaneous Rayleigh splitting occurs 22 1 . 



Figure 2. Dependence of the dimcnsionlcss chain size a on the dimensionless inverse screening length 0. The results 
are the numerical solution of Eqs. (JSJ) and {9]) for JV = 100, / 2 ^s/6 = 0.1 and r = 0.5. 



Figure 3. Variation of the dimensionless chain size a as a function of the reduced temperature r. The re- 
sults are the numerical solution of Eqs. (JSJ) and ([5]) for N = 100, / 2 ^s/6 = 0.05, and three different values of nb (a) 
first-order rod-to-globule transition via the intermediate pearl necklace structure at Kb = 0.17; (b) weak first-order 
transition without the intermediate structure at Kb = 0.35; (c) second-order continuous transition at Kb = 0.42. 
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